The study on the effect of mercury pollution on soil microorganisms around mercury mining area

In order to further explore the effects of soil mercury pollution on soil microbial diversity and community structure, soil samples were randomly collected from 2 m, 20 m, 30 m, 500 m and 650 m periphery of Wanshan mining area, as 5 different treatments. Each treatment had 4 replicates. Soil microbial DNA was extracted from 20 soil samples, and then high-throughput sequencing technology was used to analyse the structure and distribution of bacterial and fungal communities. The results showed that the number of bacterial and fungal communities in T0–T30 treatments was significantly larger than that in T500–T650 treatments at order, family and genus level. Whatever, the number of uniquely distributed bacterial and fungal communities among 4 replicates soil samples was quite different at order, family and genus level. The results of the effect on the microbial community structure showed that there were both the same dominant bacterial and fungal communities, and the different dominant bacterial and fungal communities at any classification level, moreover, the number of same dominant bacterial and fungal communities was larger than that of different dominant bacterial and fungal communities. The results of relationship between soil environment factors and bacterial and fungal community structure showed that distance (Hg2+), EC and pH had a high correlation with community structure, especially the distance factor, that is, the content of mercury in soil had the highest effects on community structure. The internal heterogeneity of soil caused significant differences in bacterial and fungal community structure, and the emergence of dominant bacterial and fungal communities was a manifestation of better adaptability to long-term mercury stress and other stresses in soil, which will provide a scientific reference for further exploring the mechanism of mercury enrichment between microorganisms and plants.


Sample collection
Five sampling points were selected at different distances around the abandoned area of Wanshan Mercury Mine in Tongren, Guizhou Province, China.Soil samples were collected at 16:44 p.m. on May 9, 2022, in summer.The T0 treatment soil samples, numbered MSO1, MSO2, MSO3 and MSO4, were collected at different points about 2 m from the sewage outlet of the mercury mining area.The T20 treatment soil samples, numbered CMSO21, CMSO22, CMSO23, and CMSO24 were collected from the wasteland of perennial ryegrass about 20 m away from the mining area.The T30 treatment soil samples, numbered GMSO301, GMSO302, GMSO303, and GMSO304 were collected from the wasteland of perennial ryegrass about 30 m away from the mining area.The T500 treatment soil samples, numbered GMSO500, GMSO501, GMSO502, and GMSO503 were collected from new cultivated land at different points about 500 m away from the mercury mining area.The T650 treatment soil samples, numbered GMSO651, GMSO652, GMSO653 and GMSO654 were collected from long-term cultivated land at different points about 650 m away from the mercury mining area, each treatment had 4 replicates.

The basic physical and chemical properties of the tested soils
The basic physical and chemical properties of soil samples were shown in Table 1.The content of basic nutrients in the soil samples was determined according to the standard analysis method 17 .In all 20 soil samples, pH was from 6.5 to 8.5; the content of available N was from 7.94 to 9.11 mg/kg, the mean of available N was T0 > T650 > T30 > T500 > T20; the content of available P was from 10.13 to 12.14 mg/kg, the mean of available P was T650 > T30 > T0 > T500 > T20; the content of available K was from 27.24 to 30.14 mg/kg, the mean of available P was T30 > T0 > T500 > T650 > T20; the content of Hg 2+ was from 45.42 mg/kg to 167.14 mg/kg, the mean of Hg 2+ was T20 > T0 > T650 > T30 > T500, Hg 2+ content of T0 treatment and T20 treatment was higher than that of T30, T500, and T650 treatment (P < 0.05), Hg 2+ content of T650 treatment was higher than that of T30 and T500 treatment (P < 0.05).

DNA extraction
Total community genomic DNA extraction was performed using a E.Z.N.A™ Mag-Bind Soil DNA Kit (Omega,M5635-02,USA), following the manufacturer's instructions.The concentration of the DNA was measured by using a Qubit 4.0 (Thermo,USA) to ensure that adequate amounts of high-quality genomic DNA had been extracted.

ITS rRNA gene amplification by PCR
Our target was the ITS1-ITS2 hypervariable region of the fungal ITS rRNA gene.PCR was started immediately after the DNA was extracted.The ITS rRNA ITS1-ITS2 amplicon was amplified using 2 × Hieff ® Robust PCR Master Mix (Yeasen, 10105ES03, China).Two universal fungal ITS rRNA gene amplicon PCR primers (PAGE purified) were used: the amplicon PCR forward primer (CTT GGT CAT TTA GAG GAA GTAA) and amplicon PCR reverse primer (GCT GCG TTC TTC ATC GAT GC).The reaction was set up as follows: microbial DNA (

16S, ITS gene library preparation, quantification, and sequencing
Hieff NGSTM DNA Selection Beads (Yeasen, 10105ES03, China) were used to purify the free primers and primer dimer species in the amplicon product.Samples were delivered to Sangon BioTech (shanghai) for library preparation using universal lllumina adaptor and index.Before sequencing, the DNA concentration of each PCR product was determined using a Qubit ® 4.0 Green double-stranded DNA assay and it was quality controlled using a bioanalyzer (Agilent 2100, USA).Depending on coverage needs, all libraries can be pooled for one run.The amplicons from each reaction mixture were pooled in equimolar ratios based on their concentration.Sequencing was performed using the lllumina MiSeq system (lllumina MiSeq, USA), according to the manufacturer's instructions.

Sequence processing, OTU clustering, Representative tags alignment and Biological classification
After sequencing, The two short lllumina readings were assembled by PEAR software (version 0.9.8) according to the overlap and fastq files were processed to generate individual fasta and qual files, which could then be analyzed by standard methods.The effective tags were clustered into operational taxonomic units (OTUs) of ≥ 97% similarity using Usearch software (version 11.0.667).Chimeric sequences and singleton OTUs (with only one read) were removed, after which the remaining sequences were sorted into each sample based on the OTUs.The tag sequence with the highest abundance was selected as a representative sequence within each cluster.Bacterial and fungal OTU representative sequences were classified taxonomically by blasting against the RDP Database and UNITE fungal ITS Database, respectively.

The distribution of bacterial and fungal communities at the order level
At the order level, the distribution of bacterial communities in the 12 soil samples of the proximal group (CR) was shown in Fig. 2.Among them, the number of co-distributed bacterial communities in 12 soil samples was 72, and the number of uniquely distributed bacterial communities in the 12 soil samples was very different.Among them, the number of uniquely distributed bacterial communities in MSO1 was 5, only 1 in MSO3 and CMSO23 soil samples, no uniquely distributed bacterial communities in the other 9 soil samples.The bacterial community distribution of 8 soil samples in the distal group (FBR) was shown in Fig. 3.The number of co-distributed bacterial communities in the 8 soil samples was 56, of which 10 uniquely distributed bacterial communities in GMSO651 soil samples, 3 uniquely distributed bacterial communities in GMSO503 soil samples, 2 in GMSO653 soil samples, and 1 in GMSO500 soil samples.The results of bacterial community distribution of all 20 soil samples showed that the number of bacterial communities in each soil sample of 12 soil samples in the proximal group (CR) was about 72, and in each soil sample of 8 soil samples in the distal treatment group (FBR) was about 58.The number of bacterial communities in T0, T20, and T30 treatments was larger than that in T500 and T650 treatments.At the order level, the distribution of fungal communities in 12 soil samples of the proximal group (CR) was shown in Fig. 4. The total number of fungal communities in all 12 soil samples was 30, but the number of uniquely distributed fungal communities to each of the 12 soil samples was very different.Among them, the number of uniquely distributed fungal communities in MSO3 soil sample was 3, 2 uniquely distributed fungal www.nature.com/scientificreports/communities in CMSO21, CMSO23, CMSO24, GMSO302, and GMSO303 soil samples, 1 uniquely distributed fungal communities in MSO2 and MSO4 soil samples, and 0 uniquely distributed fungal communities in MSO1, CMSO22, GMSO301, and GMSO304 soil samples.The distribution of fungal communities in 8 soil samples of the distal group (FBR) was shown in Fig. 5.The number of co-distributed fungi in 8 soil samples was 19, of which 9 uniquely distributed fungal communities in GMSO501 soil samples, 5 uniquely distributed fungal communities in GMSO503 soil sample, 3 uniquely distributed fungal communities in GMSO654 soil sample, 1 uniquely distributed fungal communities in GMSO500 and GMSO651 soil samples, and 0 uniquely distributed fungal communities in GMSO652 and GMSO653 soil samples.From the analysis of the number of fungal communities in all 20 soil samples, the number of fungal communities per sample in the 12 soil samples was about 31, and the number of fungal communities per soil sample in the 8 soil samples was about 22.The number of fungal communities in T0, T20, and T30 treatments was larger than that in T500 and T650 treatments.At the order level, the total number of bacterial communities in both the proximal group and the distal group was larger than that of fungal communities.

The distribution of soil bacterial and fungal communities at family level
At the family level, the distribution of bacterial communities in 12 soil samples in the proximal group (CR) was shown in Fig. 6.The total number of bacterial communities in 12 soil samples was 114, and the number of uniquely distributed bacterial communities to each of the 12 soil samples was quite different.Among them, the number of uniquely distributed bacterial communities in MSO1 soil sample was 10. 1 uniquely distributed bacterial communities in MSO2, MSO3, MSO4, CMSO21, CMSO22, CMSO23 and CMSO24 soil samples, while 0 uniquely distributed bacterial communities in GMSO301, GMSO302, GMSO303 and GMSO304 soil samples.The distribution of bacterial communities in 8 soil samples of the distal group (FBR) was shown in Fig. 7.There www.nature.com/scientificreports/were 89 bacterial communities in 8 soil samples, and the number of uniquely distributed bacterial communities in 8 soil samples was quite different.Among them, the number of uniquely distributed bacterial communities in GMSO651 soil sample was 22, 7 uniquely distributed bacterial communities in GMSO653 soil sample, 3 uniquely distributed bacterial communities in GMSO500, GMSO503 and GMSO652 soil samples, 1 uniquely distributed bacterial communities in GMSO502 soil sample, 0 uniquely distributed bacterial communities in GMSO501 and GMSO654 soil samples.From the analysis of the number of bacterial communities in all 20 soil samples, the number of bacterial communities in each sample of 12 soil samples was about 115, and the number of bacterial communities in each sample of 8 soil samples was about 92.The number of bacterial communities in T0, T20, and T30 treatments was larger than that in T500 and T650 treatments.At the family level, the distribution of fungal communities in 12 soil samples of the proximal group (CR) was shown in Fig. 8.In the proximal group (CR), the number of co-distributed fungal communities in 12 soil samples was 36, and the number of uniquely  distributed fungal communities in 12 soil samples was quite different.Among them, the number of uniquely distributed fungal communities in CMSO21 soil sample was 7, 6 uniquely distributed fungal communities in MSO4 and GMSO302 soil samples, 5 uniquely distributed fungal communities in MSO3 soil sample, 4 uniquely distributed fungal communities in MSO2 and CMSO24 soil samples, 3 uniquely distributed fungal communities in GMSO303, GMSO304, and CMSO23 soil samples. 1 uniquely distributed fungal communities in CMSO22 and GMSO301 soil samples, while 0 uniquely distributed fungal communities in MSO1soil sample.The distribution of fungal communities in 8 soil samples of the distal group (FBR) was shown in Fig. 9.The number of co-distributed fungal communities in 8 soil samples was 16, and the number of uniquely distributed fungal communities among the 8 soil samples was quite different.Among them, the number of uniquely distributed fungal communities in GMSO503, GMSO501 and GMSO502 soil samples was 22, 19 and 15 respectively, while 1 uniquely distributed fungal communities in GMSO500 soil sample.Among GMSO651, GMSO652, GMSO653 www.nature.com/scientificreports/and GMSO654 soil samples, only GMSO654 soil sample had 6, GMSO651 and GMSO653 soil samples had 1 respectively, and GMSO652 soil sample had 0. From the distribution of fungal communities in all 20 soil samples, the number of fungal communities in each sample of 12 soil samples was about 39, and the number of fungal communities in each sample of 8 soil samples was about 22.The number of fungal communities distribution in T0, T20, and T30 treatments was larger than that in T500 and T650 treatments at the family level.

The distribution of bacterial and fungal communities at genus level
At the genus level, the bacterial community distribution of 12 soil samples in the proximal group (CR) was shown in Fig. 10.The number of co-distributed bacterial communities in 12 soil samples was 173, and the number of uniquely distributed bacterial communities to each of the 12 soil samples was very different.Among them, the number of uniquely distributed bacterial communities in MSO1 soil sample was the largest, reaching 29. 1, 4, and 1 uniquely distributed bacterial communities in MSO2, MSO3, and MSO4 soil samples respectively.3, 5, 3, and 2 uniquely distributed bacterial communities in CMSO21, CMSO22, CMSO23, and CMSO24 soil samples respectively.1, 1, 0, 1 uniquely distributed bacterial communities in GMSO301, GMSO302, GMSO303 and GMSO304 soil samples respectively.The distribution of bacterial communities of 8 soil samples in the distal (FBR) group was shown in Fig. 11.The number of co-distributed bacterial communities in 8 soil samples was 101, and the number of uniquely distributed bacterial communities to each soil sample varied greatly.Among them, the number of uniquely distributed bacterial communities in GMSO651 soil samples was 62, and 12, 14, 1 uniquely distributed bacterial communities in GMSO652, GMSO653, and GMSO654 soil samples respectively.4, 6, 3 and 3 uniquely distributed bacterial communities in GMSO500, GMSO501, GMSO502 and GMSO503 soil samples respectively.According to the distribution of bacterial communities in all 20 soil samples, the number of bacterial community distribution in each sample of 12 soil samples was about 177, while the number of bacterial communities in each sample of 8 soil samples was about 114, the number of bacterial communities in T0, T20, and T30 treatments was larger than that in T500 and T650 treatments.At the genus level, the distribution of fungal communities in 12 soil samples of the proximal group (CR) was shown in Fig. 12.The number of co-distributed fungal communities in 12 soil samples was 21, but there was a big difference in the number of uniquely distributed fungal communities among the 12 soil samples.The number of uniquely distributed fungal communities in MSO1, MSO2, MSO3 and MSO4 soil samples was 2, 13, 14, 21 respectively.22, 5, 12, 16 uniquely distributed fungal communities in CMSO21, CMSO22, CMSO23 and CMSO24 soil samples respectively.6, 10, 11, 12 uniquely distributed fungal communities in GMSO301, GMSO302, GMSO303 and GMSO304 soil samples respectively.The distribution of fungal communities in 8 soil samples of the distal group (FBR) was shown in Fig. 13.The total number of fungal communities in the 8 soil samples was 12, but the number of uniquely distributed fungal communities among the 8 soil samples varied greatly.Among them, the number of uniquely distributed fungal communities in GMSO500, GMSO501, GMSO502, GMSO503 soil samples was 8, 50, 44, 55, respectively, and 14, 6, 2, 11 uniquely distributed fungal communities in GMSO651, GMSO652, GMSO653, GMSO654 soil samples respectively.From the distribution of fungal communities in all 20 soil samples, the number of fungal communities each sample in 12 soil samples was about 33, and the number of fungal communities per sample in the 8 soil samples was about 35, The number of fungal communities in T0, T20, and T30 treatments was close to the number of fungal communities in T500 and T650 treatments.www.nature.com/scientificreports/

Beta diversity analysis
The beta diversity analysis of bacterial communities among 20 soil samples was shown in Fig. 14 at the OTU level.From the bacterial community structure of 20 soil samples in T0-T650 treatments, except for the MSO3,  www.nature.com/scientificreports/CMS021, GMSO502 soil samples, the bacterial community composition of the other 18 soil samples was nearer to that of the soil samples in their respective treatment.Based on the similarity of bacterial communities of soil samples among T0-T650 treatments, from the perspective of phylogenetic distance, the results of the occurrence of bacterial communities in 5 treatments of soil samples from near to far was T20, T0, T500, T30, T650.By comparing the similarity of bacterial communities between CR and FBR group soil samples, the results showed that bacterial communities of T650 treatment in FBR group soil samples had lower similarity with T0, T20, T30 treatments in CR groups soil samples, but bacterial communities of T500 treatment in FBR group soil samples had higher similarity with T0, T20, T30 treatments in CR groups soil samples.

The relationship between soil environment factors with bacterial and fungal abundance
The relationship between soil environment factors and bacterial community abundance of CR and FBR group soil samples was shown in the Fig. 22.The results showed that among 6 soil environmental factors, distance, pH and EC had very important effects on the distribution of bacterial communities, In particular, the distance had the greatest effect, indicating that mercury in the soil played a dominant role in the distribution of bacterial       communities.The bacterial community abundance of FBR group soil samples was more affected by distance than that in CR group soil samples.
The relationship between soil environment factors and fungal community abundance of CR and FBR group soil samples was shown in the Fig. 23.The results showed that among the 6 soil environmental factors, distance, EC and pH had very important effects on the distribution of fungal communities, In particular, the distance had the greatest effect, indicating that mercury in the soil played a dominant role in the distribution of fungal communities.The fungal community abundance of FBR group soil samples was more affected by distance than that of CR group soil samples.While the fungal community abundance of CR group soil samples was more affected by pH, and EC than that of FBR group soil samples.

Discussion
Mercury is one of the most common heavy metal pollutants in the environment.Soil mercury pollution is often closely related to the mining of surrounding mercury-related metal mines and the smelting of mercury-related metals.Once soil is polluted by mercury, it will be transferred to plant roots, stems, leaves and other parts through plant enrichment.When the soil mercury content continues to rise, it will seriously affect the normal growth and development of plants 18,19 .Meanwhile, it will have adverse effects on the soil ecosystem, especially on soil  www.nature.com/scientificreports/microorganisms and soil animals 20 .Previous studies have confirmed that the concentration of Cd, Zn, Pb and Cu in soil had a strong correlation with soil microbial biomass 21 , and different microbial groups had different tolerance to different heavy metals 22 .The reason why heavy metals had adverse effects on microbial diversity could be that high concentration heavy metals had a certain restrictive effect on microbial cell metabolism and other functions, which ultimately led to a decrease in microbial diversity 23 .In this work, the results of bacterial and fungal community composition in soils with different mercury contents at different distances from mercury mining areas showed that the number of bacterial and fungal communities in 20 soil samples increased significantly from the order level to the genus level.The number of uniquely distributed bacterial and fungal communities to each soil sample of the 20 soil samples showed a significant increase trend from the order level to the genus level.This result showed that even in a certain mercury stress environment, the number of bacterial and fungal communities in the soil increased with the decline of the grade level, which was not significantly affected by soil mercury and other factors.Meanwhile, the alpha diversity of bacterial and fungal communities were analyzed by using reads, OTUs, Shannon index and Simpson index (showed in Tables 2 and 3), our results indicated that there were large differences in the diversity of bacterial and fungal communities among 4 replicates soil samples, from T0 to T650 treatments, the diversity of bacterial communities decreased gradually, however, the diversity of fungal communities increased gradually.This results indicated that in the outdoor natural environment, the diversity distribution of bacterial and fungal communities at different points within the same range was affected by the internal factors of the soil.The diversity of bacterial and fungal communities was affected by many factors, resulting in the difference in the number of uniquely distributed bacterial and fungal communities to individual soil samples 24 .It also indicated that there was a certain heterogeneity in the soil around the mercury mining area.Therefore, it can be speculated that mercury press will lead to a decline in soil microbial diversity to a certain extent, but the relationship between the two was not a simple linear relationship, also affected by other factors within the soil.Meanwhile, the results on beta diversity analysis of bacterial and fungal communities showed that there were some differences between bacterial and fungal communities, the similarity of bacterial communities from far to near was T20, T0, T500, T30, and T650 treatment soil samples, the similarity of fungal communities from near to far was T0, T30, T20, T650, and T500 treatment soil samples.The similarity of bacterial communities and fungal communities in different soil samples was significantly different, and the main reason for this result was related to the internal heterogeneity of soil.
The results of this study showed that the number of bacterial communities in T0 and T20 treatments was much larger than that in T30, T500, and T650 treatments at order, family and genus level.Also lying in fungal communities.That is, the number of bacteria and fungal communities decreased with decline of mercury concentration at order, family and genus level.This result was inconsistent with the research results of Shan and coworkers 25 .This could be caused by the difference of mercury concentration or experimental environment.Previous studies belonged to the artificial simulation laboratory soil environment.This study was the natural soil environment of mercury mining area, various bacterial and fungal communities reached an equilibrium state after a long period of comprehensive influence of various factors inside the soil.The real existence of these bacterial and fungal communities was a good adaptation to mercury-containing soil environment for a long time.The reason for this result could be related to the soils of different types.Maliszewska and coworkers 26 found that the compounds of Cu and Hg showed the strongest inhibitory effect on microbial proliferation, and the order of toxicity of HgCl 2 to microorganisms was bacteria > fungi > actinomycetes.The results showed that the inhibitory effect of mercury on microbial proliferation was fluctuated, and the inhibition of microorganisms decreased with the prolongation of culture time 26 .In this work, the soil around the mercury mining area had been subjected to high concentration mercury pollution stress for a long time.The increase in the number of bacterial and fungal communities under higher concentration mercury environment could be the best evidence that the inhibition of bacterial and fungal communities decreased with time.
The structure of bacterial, fungal and communities of actinomycetes in soil was easily affected by environmental variables around the soil, such as nutrient cycling and heavy metal concentration changed 27,28 .In this work, the effects of mercury concentration on the structure of bacterial and fungal communities at the order, family and genus levels were analyzed.The results indicated that mercury had a certain effect on the structure of bacterial and fungal communities, and the tolerance of different bacterial and fungal communities to mercury pollution was also very different, which was consistent with the results of Prasad 29 .Our results indicated that the dominant bacterial communities in T0-T650 treatments were different at order, family and genus level.At the order level, Sphingomonadales and Sphingobacteriales were the dominant bacterial communities, the relative abundance remained at about 10%; at the family level, for Sphingomonadaceae and Chitinophagaceae; at the genus level, for Gp6(norank_Acidobacteria_Gp6), Sphingomonas,Gp4 (norank_Acidobacteria_Gp4).Similar results were also found in fungal communities, at the order level, Agaricales and Sebacinales were the dominant fungal communities; at the family level, for Sebacinaceae and Hygrophoraceae; at the genus level, for Hygrocybe and Sebacina.Regardless of the classification level, bacterial and fungal communities with higher relative abundance were the performance of long-term adaptation to higher mercury soil environment, indicating that they had strong tolerance to heavy mercury stress.This was similar to the conclusion of Feris 24 and Gillan et al. 30 .Meantime, at order, family and genus level, the dominant bacterial and fungal communities in different positions of mercury mining area were quite different, which was closely related to the heterogeneity of many factors in the internal environment of soil.The relationship between soil environment factors and bacterial and fungal community abundance was analyzed by using db-RDA method at OTU level, the results indicated that distance(Hg 2+ concentration), EC, and pH three soil factors had greater effects on the relative abundance of bacterial and fungal communities than available N, P, K.Among distance(Hg 2+ concentration), EC, and pH, distance(Hg 2+ concentration) had the greatest effects on the relative abundance of bacterial and fungal communities.This results fully demonstrated that these dominant bacterial and fungal communities had much stronger resistance to mercury pollution stress.

Figure 1 .
Figure 1.Location map of study area.

Figure 2 .Figure 3 .
Figure 2. The distribution of bacterial communities in 12 soil samples of the proximal group (CR) at the order level.

30 Figure 4 . 19 Figure 5 .
Figure 4.The distribution of fungal communities in 12 soil samples of the proximal group (CR) at the order level.

Figure 6 .
Figure 6.The distribution of bacterial communities in 12 soil samples of the proximal group (CR) at the family level.

Figure 7 .
Figure 7.The distribution of bacterial communities in 8 soil samples of the distal group (FBR) at the family level.

Figure 8 . 15 Figure 9 .
Figure 8.The distribution of fungal communities in 12 soil samples of the proximal group (CR) at the family level.

Figure 10 .
Figure 10.The distribution of bacterial communities in 12 soil samples of the proximal group (CR) at the genus level.

Figure 11 . 2 M 21 Figure 12 .
Figure 11.The distribution of bacterial communities in 8 soil samples of the distal group (FBR) at the genus level.

12 Figure 13 .
Figure 13.The distribution of fungal communities in 8 soil samples of the distal group (FBR) at the genus level.

Figure 17 .
Figure17.The relative abundance of fungal communities in 20 soil samples at the order level.

Figure 18 .
Figure 18.The relative abundance of bacterial communities in 20 soil samples at the family level.

Figure 19 .
Figure 19.The relative abundance of fungal communities in 20 soil samples at the family level.

Figure 20 .
Figure 20.The relative abundance of bacterial communities in 20 soil samples at the genus level.

Figure 21 .
Figure 21.The relative abundance of fungal communities in 20 soil samples at the genus level.

Figure 22 .
Figure 22.Redundancy analysis (db-RDA) between soil environment factors and bacterial abundance of CR and FBR group at OTU level.

Figure 23 .
Figure 23.Redundancy analysis (db-RDA) between soil environment factors and fungal abundance of CR and FBR group at OTU level.

Table 1 .
Physiochemical properties of the tested soils.Data are means土standard deviation.The different case letters indicate that the means are significantly different among soils (P < 0.05) with Duncan test.

Table 2 .
Statistical table of bacterial communities alpha diversity.

Table 3 .
Statistical table of fungal alpha diversity.